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ABSTRACT 

We study a cubic lattice gas model for nuclear matter where each lattice site can be either 
occupied, by one proton or one neutron, or unoccupied. A nearest-neighbor interaction of the 
form — J2<ij> Jij T ziT z j is assumed. Our model is an isospin-1 Ising model, with r z = (1,0,-1) 
representing respectively (proton, vacancy, neutron). A kinetic-energy term has been included 
in our model. Under the Bragg- Williams mean field approximation our model exhibits the 
existence of a dense phase (liquid- like) and a rare phase (gas- like). The nuclear-matter p-v 
isotherms given by our model are discussed. 
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It is believed that nuclear matter should exhibit a quark deconfinement phase tran- 
sition at high temperatures {k B T ~ 200 MeV) and/or high densities (p ~ 10po), p 
being the saturation density. At low temperatures (&bT ~ 15 — 20 MeV) and low 
densities (p ~ po) a liquid-gas phase transition of nuclear matter is expected to take 
place [|T|]. The underlying theoretical frameworks treating these two types of phase 
transitions have so far been very different. 

For the former one generally uses lattice gauge QCD and Monte Carlo simulation 
[0. The space-time structure here is a lattice. But for the latter, it has been a long 
tradition to employ standard many-body theories such as Hartree-Fock approximation 
||, |J, real time Green's function method H ||, and ring-diagram summation method 
0. The space-time structure here is a continuous manifold. There appears to be a 
space-time disparity and one would like to ask the following questions: Is it necessary 
for nuclear systems to have so very different space-time geometry structures, namely 
a lattice and a continuum, in different energy and density domains? Can we set our 
theory on a unified footing so that a common lattice space-time manifold may be used 
to treat not only phase transitions in high-temperature and high-density regions but 
also those with low temperatures and low densities? 

The above concern has motivated us to explore new schemes for studying nuclear 
matter in the low-temperature and -density regions. About forty years ago Lee and 
Yang || suggested a model of lattice gas where gas atoms are seated on a lattice. They 
mapped the problem of a lattice gas with one type of atom into an Ising model for spin 
half particles. They succeeded in describing a liquid-gas phase transition for atomic 
systems. We would like to generalize this model to nuclear matter which consists of 
two types of nucleons: proton and neutron. Here we consider a lattice with each site 
either being vacant or occupied by a proton or a neutron. One is tempted to associate 
Ising spin 1 to a proton, -1 to a neutron, and to a vacancy. A preliminary account 
of our approach has been reported 0. 

Nuclear forces may be understood based on one boson exchanges . The intermediate 
bosons include isoscalar bosons (u,a,...) and isovector bosons (7r,p, ...). In this way, 
effective nucleon-nucleon(NN) interactions may be written in terms of a number of 
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standard components such as the spin-spin term a ■ a , isospin-isospin term r • r, and 
the spin-isospin term a ■ a r ■ r. If we average over the spin and spatial variables, 
we would obtain an effective interaction which depends on isospin only. The radial 
part of the NN interaction, to a reasonable approximation, can be represented by a 
nearest-neighbor square well potential with a repulsive hard core. Thus we are led to 
conjecture the following interaction hamiltonian for nuclear matter: 

Hint ^ t Jij^zi^zj (1) 

<ij> 



where 



1 for proton 

i~ z i = { for vacancy (2) 
-1 for neutron. 



Note the above is a nearest-neighbor interaction, as indicated by the summation index 
< ij >, namely interaction exists only between adjacent neighbors. Jij are interaction 
strength parameters. In fact we allow for two such parameters only; J s for pp and nn 



interaction pairs and Jd for pn pairs. Clearly the above is a spin 1 Ising model [jTO], [IT 
We rewrite eq.(l) as 

H mt = -J S (N ++ + AL_) + J d N + _ (3) 

where N ++ ,N represent respectively the nearest neighbor pairs of proton-proton, 

neutron-neutron, etc. Note that in our model there is no interaction between vacancy 
and nucleon. 

Now we proceed to calculate the thermodynamic quantities of our system. To 
facilitate the calculation, we introduce the following variables. Let N denote the total 
number of lattice sites, and N + , N_, N denote respectively the number of proton, 
neutron and vacancy sites. We introduce 

N n N.—N 
R= T> S= N "" ' ^ = ^+ + ^- + ^0 (4) 
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where R and S represent the relative emptiness and proton-neutron asymmetry re- 
spectively. Note that in our model the nuclear-matter density p is proportional to 
(1-R). 

The spin- 1/2 lattice-gas models for atomic systems have been used with remarkable 
success in describing phase transitions. But they seem to have not, to our knowledge, 
been applied to calculate the pressure-volume isotherms. It has been suggested 
that in order to describe these isotherms one needs to add the ideal-gas pressure to the 
lattice gas grand potential. The ideal gas pressure comes from the kinetic energy. Hence 
it would seem to be preferrable to include the the kinetic energy in the hamiltonian, 
from the beginning, rather than adding the ideal-gas pressure to the grand potential. 
The above consideration has motivated us to include a kinetic-energy term in our 
model, namely we employ a hamiltonian for the lattice gas of the form 

H gas = -J S (N++ + AL_) + J d N + „ + Nk{\ - R) 5/3 (5) 

where, guided by the Fermi gas model, we have assumed the kinetic energy per particle 
to be proportional to p 2//3 . k is a constant which we shall discuss later. The introduction 
of the kinetic-energy term is very important for our model, as we shall see soon. 

As an initial investigation let us adopt the Bragg- Williams mean field approximation JT2 
namely 

N 2 + N ++ Nl iV__ iV 2 _ N 00 



N 2 JV7/2' N 2 JV7/2' N 2 N 1 /2 
where 7 denotes the number of nearest neighbors of any given site, and N'~f/2 is the 
total number of pairs. For three dimensional simple cubic lattice, 7 = 6. From the 
constraints -fN + = 2N ++ + iV + _ + iV +0 , 7AL = 2AL_ + iV + _ + AL , 7^ = 2A^ 00 + 
iV +0 + iV_ and iV = iV + + iV_ + N , we obtain 

N + _ = 7 (iV + + iV_) - (N ++ + iV__) + iVoo - jN/2. (7) 

We can now rewrite our hamiltonian of eq.(5) as 

H gas (R, S, N) = -dNS 2 - C 2 N{1 - Rf + Nk(1 - R) 5/3 (8) 
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where 



n {J s + Jdh n {J* - Jdh / q n 

Cl 4 ' 62 4 • ^ 



The grand partition function of our system is 

Qg = T,9(R, S,N)exp(-H 9as /k B T), 

R,S 

H gas = H gas (R, S, N) - hNS - XN(1 - R) (10) 
The multiplicity factor |[3] g is 

3 S ' N ^ = N \N + \NJ. = (NR)\[N(l-R + S)/2]\[N(l-R-S)/2]\' ^ 

In H gas we have included two Lagrange multipliers, h and A. In our grand partition 
function, we sum over all possible R and S values. Thus we have neither a definite num- 
ber of nucleons nor a definite proton-neutron assymmetry. The role of these Lagrange 
multipliers is to control the average values of R, denoted as R, and of S, denoted as S. 
For instance for symmetric nuclear matter we need >S=0, and this may be attained by 
varying h. Similarly, different R values may be obtained by varying A. 

Since our final goal is to study the thermodynamic limit (N — > oo), we can replace 
the sum in the partition function Qg by its most dominant term JO?] (assuming the 
dominant term to be non-degenerate). Using Stirling's formula one obtains 



k B TlnQ G = —CiNS 2 — C 2 N(1 — R) 2 — hNS 



5/3 



-XN(l -R) + kN(1 - R) 
+Nk B T[R \nR+ ^~ R + S ^ ln (! _ R + S ) 



+ (1 — ^_^li n (i_ jR _,s)_(i_ jR )i n 2] (12) 
2 

with the subsidiary extremum conditions 

d(-k B T\nQ G ) Nk B T 1-R + S lorqj _ h]Af n h ^ 

ds = — ln T^R^s- [2ClS + h]N = ' (13) 



Kuo, Ray, Shamanna and Su / Isospin Lattice Gas Model and Nuclear matter 6 

d(-k B T\nQ G ) _ Nk B T 4R* 

OR " "^~ ln (1-R)*-S* +2NC2{1 - R) 

- ^kN(1-R) 2/3 + X = 0. (14) 

It is seen that h — and S = is a special solution which represents a symmetric 
nuclear matter. In this case eq.(13) becomes an identity and eq.(14) reduces to 

or c; 

k B T\n - + 2C 2 (1 -R)- -k(1 - i?) 2 / 3 + A = 0. (15) 

1 — rt o 

The analysis of this equation gives us the most important results of our present 
paper. It shows that below a certain temperature T c we have two phases, one is dense 
(liquid-like) and the other is rare (gas-like). Above this temperature there is only one 
phase. The existence of a liquid-gas phase transition together with the determination of 
its phase diagram on the basis of a simple model that assumes only a phenomenological 
two body, nearest-neighbor interaction is quite remarkable. It seems to call for serious 
attention and further study. 

Let us consider the A=0 case first. In this case we rewrite eq.(15) as 

X (R,T) = f(R)-g(R,T) = 0, 

We note that f(R) is a monotonically increasing, unbounded (at R — > 0+ , and 
R — > 1— ) function of R having one point of inflection at R = 1/2. However g(R,T) 
is a bounded function in the same domain. Hence g(R,T) must intersect f(R) at least 
once, i.e., eq.(16) must have at least one solution. In addition g(R,T) has a negative 
curvature and one maximum for R G [0, 1]. Thus there is a possibility of having more 
than one point of intersection with f(R) below a suitable temperature T c . 

In Fig. 1 we display some typical behaviours of f(R) and g(R,T). As shown, f(R), 
denoted by the solid line, approaches to — oo at R=0 and to oo at R=l. g(R,T) is 
a concave-downward curve. f(R) intersects the R axis at R — 1/3. g(R,T) intersects 
the R axis at two points: One is at R=l which is independent of the values of the 
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parameters C 2 , k, and ksT, and the second point of intersection depends on the ratio 
a = 5/t/6C 2 and importantly is independent of temperature T. 

For low temperatures, the curves f(R) and g(R,T) may have three intersection 
points, as denoted by A,B, and C for the T=8 case. (Note that we use the convention 
of ks = 1.) It is readily checked that the middle intersection point, i.e. B, corresponds 
to a minimum of lnQ G an d hence it is not a physical solution. The intersection points 
A and C are the physical solutions. Note that g(R,T) is inversely proportional to the 
temperature. Hence as T increases, the right side of it sweeps down while its left side 
sweeps up, as shown in the figure. At some critical temperature T c the curves f(R) 
and g(R,T c ) become tangent to each other. And afterwards the curves f(R) and g(R,T) 
have only one intersection point. g(R,T) eventually becomes a horizontal line as T 
approaches 00, as indicated by the T = 00 dotted line in the figure, and it obviously 
has just one intersection point with f(R). 

In Fig. 2 we plot the solutions of eq.(16). As seen, we have two solutions, as 
indicated by points A and C of Fig. 1, for T < T c , and for T > T c we have only one 
solution. This figure suggests that below T c we have three regions, the dense(liquid- 
like) phase, the rare(gas-like) phase and the coexistence phase in between. The physical 
meaning of the above results may become clearer by examining the pressure-volume 
isotherms. 

Before doing this, we should probably discuss the parameters C 2 amd k which have 
entered into our calculations. To have an attractive nearest-neighbor interaction, we 
have J s > and J d < 0. In this case C 2 is positive. The magnitude of J s and J d is 
probably comparable to the average potential energy in nuclear matter, which is about 
-40 MeV. Hence we have picked C 2 = 125 MeV, recalling that C 2 has been defined 
in eq.(9). The parameter k may be estimated from the average kinetic energy given 
by the Fermi gas model. Assuming a lattice spacing of 1.5 fm and taking the nucleon 
mass as 940 MeV, k is obtained as about 125 MeV. Hence as a general guideline, the 
parameters C 2 and k should both be not too far from 125. The results of Figs.l and 
2, and 3 to be presented later, are all obtained with C 2 =125 and K=1.2aC 2 with a 
specially chosen a as discussed below. 
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Near the critical point, there is a subtle dependence of the solutions of eq.(16) on the 
ratio a = jj^k As indicated by points A, B and C of Fig. 1, we have three intersection 
points below T c . It is readily checked that by choosing a = (2/3) 1 / 3 , these three 
intersection points all merge together at the critical point. Then the phase diagram 
near T c has the smooth shape as shown in Fig. 2. If one uses a slightly different ratio, 
the merging would generally take place in two steps, first involving two intersection 
points and then the third. This will lead to a phase diagram with a "cusp" shape near 
the critical point, which may seem to be rather unconventional. By intuition we feel 
that it is more reasonable to have a smooth phase boundary, and hence we have chosen 
the above a. With this a and C 2 =125, we have ^=131.037, which has been used in 
obtaining the results presented in our Figs. 1 to 3. 

For fixed T, the solutions of eq.(16) determines the R values where the grand 
partition function has a maximum contribution. In fact this maximum term is an 
overwhelmingly dominant one. (For finite lattices such as that with N=10 6 , numerical 
simulation has shown that the magnitude of the maximum term is typically ~ 10 100 
while for all the other terms it is ~ 10 20 .) Hence the R values given by eq.(16) are just 
the average values of R, as denoted previously by R, for the system at temperature 
T and A=0. The R values for A ^ are given by the solution of eq.(15). Then the 
pressure at various densities and temperatures is given by 

p(R,T) = ^f In Q G 

op 1 p 

= C 2 (l - R) 2 - k(1 - Rf 3 - k B TR\n(—^) - k B T ln( — -) + A(l - R). (17) 

Let us define the specific volume v as (1 — Ry 1 . Then we can calculate the p-v 
isotherms using eqs.(16) and (18). Our results are shown in Fig. 3. 

Our figure seems to contain some interesting features. For a given temperature T 
the isotherm is obtained as a parametric plot of specific volume versus pressure with A 
its generating parameter. For T >T C one gets a single smooth curve by varying A. The 
boundary, drawn as a solid line at v~ 1.5, corresponds to A = 0. The isotherm to the 
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left of this boundary is obtained with positive A, while to the right of this boundary the 
isotherm is obtained with negative A. The critical temperature obtained is T c = 18.4. 

For T below T c , it is of interest that no isotherms are obtained in the intermediate 
region, i.e. region I. For A < 0, we get a gas-like isotherm starting from the high specific 
volume tail of the phase boundary. And for A > we get a liquid-like isotherm (very 
high compression modulus) starting from the lower specific volume edge of the phase 
boundary. The above is because when T < T c , eq.(15) has no solutions in region I, the 
coexistence region. With the introduction of an infinitesimal A the system chooses one 
of the two values of R admissible for the given T < T c , depending on the sign of A. This 
phenomena is reminisent of the spontaneous symmetry breaking in ferromagnetism. 

There is another point which may be mentioned. As indicated in the figure, for 
a given isotherm the pressure at the liquid boundary and that at the gas boundary 
appear to be equivalent to each other, such as the apparent equivalence between the 
pressures of the T=12 isotherm at the two boundaries. We have examined this type 
of apparent equivalence, and have found that these two pressures are "exactly" equal 
to each other, within the accuracy of our computer. The structure of our calculated 
isotherms seems to strongly support that for T < T c our model gives a liquid phase, a 
gas phase and a coexistence phase with the boundary indicated by the solid line. And 
for T > T c the distinction between the liquid and gas phases disappears. 

As a conclusion, let us state the following points. We have studied a simple lattice 
gas model for nuclear matter, where each lattice site can be either vacant or occupied, 
by one proton or by one neutron. A hamiltonian consisted of a nearest-neighbor inter- 
action and a kinetic-energy term is assumed. The partition function is then calculated 
with the Bragg- William approximation. Some rather encouraging results have been 
obtained. A phase diagram consisting of liquid, gas and coexistance phases is obtained 
from our model. And the p-v isotherms given by our model are surprisingly similar 
to those given by the van der Waals theory, except for the difference that for T < T c 
our isotherms do not have the metastable states in the coexistence region. Hence with 
our model one does not need to determine the phase boundary by way of a Maxwell 
construction. The liquid-gas critical temperature obtained by us is T c = 18.4 for sym- 
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metric nuclear matter, which is fairly close to the results given by earlier calculations 



Comparing with earlier lattice gas models []12| , a new ingredient of the present 
model is the introduction of the kinetic energy term. We recall that to have phase 
transitions we need the g(R,T) curve to have three intersection points with f(R) for 
T < T c . Whether this happens or not clearly depends on a (= 5^/66*2). There is a 
wide range of values of a, for which this could happen. But if we don't have the kinetic 
energy term, then k = and a = 0. And in this case the above requirement can not 
be met, and our model would have no phase transitions. Hence the inclusion of the 



kinetic energy term is important for our model. In fact in a previous paper [14]] we 
have studied a lattice-gas model for nuclear matter, using the same lattice Hamiltonian 
as Eq. (5) but without the inclusion of the kinetic energy term; the results were clearly 
unsuccessful. 

We have adopted a major approximation, the Braggs- Williams approximation, in 
the present work. The accuracy of this approximation remains to be investigated. More 
accurate calculations may be performed, using for example the Bethe-Peierls approxi- 
mation [EJ, or Monte Carlo simulations which have been extensively in lattice-gauge 
and Ising-model calculations. With the Bragg- Williams approximation, the calculated 
phase boundary near the critical point can be either a smooth shape or a cusp shape. 
It should be of interest to see what would be the phase-boundary shape given by such 
more advanced methods. 
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FIGURE CAPTIONS 
Fig.l Graphical solution of eq.(16). 

Fig. 2 Nuclear matter phase diagram given by our model. 
Fig. 3 Nuclear matter p-v isotherms. 



